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The electronic conductance of a benzene molecule connected to gold electrodes via thiol, thiolate, 
and amino anchoring groups is calculated using nonequilibrium Green functions in combination with 
the fully selfconsistent GW approximation. The calculated conductance of benzenedithiol and ben- 
zenediamine is five times lower than predicted by standard density functional theory (DFT) in very 
good agreement with experiments. In contrast, the widely studied benzenedithiolate structure is 
found to have a significantly higher conductance due to the unsaturated sulfur bonds. These findings 
suggest that more complex gold/thiolate structures where the thiolate anchors are chemically pas- 
sivated by Au adatoms are responsible for the measured conductance. Analysis of the energy level 
alignment obtained with DFT, Hartree-Fock and GW reveals the importance of self-interaction cor- 
rections (exchange) on the molecule and dynamical screening at the metal-molecule interface. The 
main effect of the GW self-energy is to renormalize the level positions, however, its influence on the 
shape of molecular resonances also affects the conductance. Non-selfconsistent Go Wo calculations, 
starting from either DFT or Hartree-Fock, yield conductance values within 50% of the selfconsistent 
GW results. 

PACS numbers: 73.63.-b,73.40.Gk,85.65.-l-h 
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I. INTRODUCTION 

The problem of first-principles calculation of elec- 
tronic conduction in molecular systems is of longstand- 
ing interest. Over the last decade, advances in exper- 
imental techniques have allowed for fundamental stud- 
ies of electron transport through few or even individ- 
ually contacted molecules. The transport mechanisms 
observed in molecular junctions range from ballistic^ 
over off-resonant tunneling^"^ to the strong correlation 
Kondo and Coulomb blockade regime^''^ to vibration as- 
sisted hopping^. The former two belong to the phase- 
coherent transport regime characteristic of relatively 
short molecules (as opposed to molecular wiresS) with 
good " chemical" contact to the electrodes and is the main 
focus of the present work. 

The lack of control over the atomistic details of the 
metal-molecule interface introduces a strong statistical 
element in measurements on single-molecule junctions 
which masks the relation between atomic structure and 
measured conductance. For example, published conduc- 
tance values for the gold/benzenedithiol model system 
vary by several orders of magnitude, although recent in- 
dependent studies of this system seem agree on a value 
around O.OlGo (Go = 2e^//i)^"-. It is, however, not 
clear which structure is responsible for this " typical" con- 
ductance and recent experimental and theoretical work 
point to complex gold/thiolate structures involving two 
molecules bonding the same Au adatom^^"^'^. Scanning 
tunneling microscope experiments in solution have shown 
that the use of amino rather than thiol anchoring groups 
leads to more well defined junction properties^ and a con- 
ductance around (or just below) O.OlGo has also been 



reported for the gold/benzenediamine junction?'- . This 
was, however, not confirmed by independent break junc- 
tion experiments in vacuunii^. 

The uncertainties related to the junction atomic struc- 
ture renders theoretical benchmark calculations more im- 
portant and more challenging at the same time. Here 
progress has been hampered by the inability of conven- 
tional density functional theory (DFT) methods, which 
for long have been the workhorse and the state of 
the art for quantum transport calculation a^^'^^ , to re- 
produce the conductance measured for even the sim- 
plest molecular tunnel junction a^^'^" . As a consequence, 
the most succesful studies have focused on qualita- 
tive trends in the dependence of conductance on e.g. 
molecular conformational, molecular lengtb^^, side group 
functionalizationsS^, or have focused on properties inde- 
pendent of the numerical value of the conductance like 
molecular vibrationa^ii^. The main shortcoming of the 
DFT approach has been attributed to its band gap prob- 
lem, i.e. the fact that DFT tends to underestimate energy 
gaps^, and therefore overestimates the conductance. At- 
tempts to overcome this problem within a single-particle 
framework have mainly focused on self-interaction cor- 
rection schemes^l"— . 

The well known success of the many-body GW 
method^ for quasiparticle (QP) bandstructure calcu- 
lations has recently inspired its application to (simpli- 
fied) transport problems^ir— . The fact that the GW 
approximation succesfully describes systems with highly 
diverse screening properties ranging from metals^ over 
semiconductors^"-^ to molecules--'--^ is essential for a 
correct description of metal-molecule interfaces where 
the electronic character changes from metallic to in- 
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sulating over a few angstrom. In particular, screen- 
ing by the metal electrons can have large influence on 
the QP energies of the adsorbed molecule^"— - an ef- 
fect completely missed by both local and hybrid den- 
sity functionals^. This has recently motivated the 
use of semi-empirical schemes for correcting the DFT 
eigenvalues by a scissors operator prior to transport 
calculationa^^i^iS. While such schemes can be justified 
for weakly coupled molecules, they become uncontrolled 
in the relevant regime of covalently bonded molecules 
where the screening effects mix with charge transfer and 
hybridizatioit^. 

In this work we combine nonequilibrium Green func- 
tion methods for electron transport with the fully selfcon- 
sistent GW approximation for exchange and correlation 
to establish a theoretical benchmark for the electronic 
structure and conductance of gold/l,4-benzenedithiolate 
(BDT), -benzenedithiol (BDT-(-H) and -benzenediamine 
(BDA) molecular junctions. We find a conductance of 
0.0042Go for BDA and O.OlOGo for BDT+H in very 
good agreement with experimental data. In compari- 
son, the conductance obtained from DFT is about five 
times higher while non-selfconsistent GqWo calculations 
produce conductances within 50% of the selfconsistent 
GW result. We argue that the BDT-I-H structure can 
be viewed as a simple model of recently proposed RS- 
Au(I)-SR gold/thiolate structures involing two molecules 
attached to the same Au adatom. The conductance of a 
simple BDT molecule between Au(lll) surfaces is pre- 
dicted to be on the order of IGq by both DFT, Hartree- 
Fock and GW. The origin of the high conductance is 
due to an unsaturated sulfur p orbital with energy just 
below the Fermi energy. In the BDT+H and RS-Au(I)- 
SR structures, the sulfurs are fully passivated and the 
p-orbital moves away from Ep leading to an effective de- 
coupling of the C6II4 moiety from the gold electrodes. 
We find that the main effect of the GW self-energy is to 
shift the molecular levels and can be modelled by a sim- 
ple scissors operator. However, the energy dependence 
of the GW self-energy may also affect the shape of the 
transmission resonances and this can change the conduc- 
tance by almost a factor of two. 

Most implementations of the GW method invoke one 
or several technical approximations like the plasmon pole 
approximation, neglect of off-diagonal matrix elements 
in the GW self-energy, analytic continuations from the 
imaginary to the real frequency axis, neglect of core 
states contributions to the self-energy, neglect of self- 
consistency. The range of validity of these approxima- 
tions has been explored for solid state systems by a num- 
ber of author a^^'^"^" — , however, much less is known about 
their applicability to molecular and metal-molecule sys- 
tems. For this reason our implementation of the GW 
method avoids all of these approximations and as such 
represents an exact treatment of the GW self-energy 
within the space of the employed atomic orbital basis 
set. 

Although we compare our results to experimental data 



and discuss them in relation to the possible atomic struc- 
ture of the junctions, we stress that the main purpose of 
this study is the benchmarking of quantum transport cal- 
culations for specific, idealized junction with particular 
focus on the role of electronic correlation effects. 



II. METHOD 

We consider a quantum conductor consisting of a 
molecule connected to left (L) and right (R) electrodes. 
We shall assume that outside a certain region containing 
the molecule and part of the electrodes (the "extended 
molecule" ) , the electron-electron interactions can be de- 
scribed by a mean field potential. The current through 
the molecule is then given by^i^ 

^ = T^y[{rL-rii)G<+{fLrL-fRrR){G--G-)]dE 

(1) 

where the energy dependence of all quantities has been 
suppressed for simplicity. In this equation G is the Green 
function matrix of the extended molecule evaluated in 
a localized basis, F^/^ — ^[S^/^ — ^l/b) 
pling strength between the extended molecule and the 
left /right electrode, and fL/R are the Fermi Dirac dis- 
tribution functions of the two leads. Our implemen- 
tation applies to the general case of a finite bias volt- 
age, but in this work we focus on the zero bias conduc- 
tance which can be expressed in terms of the transmission 
functioit^i^ 

T{E) = Tr[G'-{E)rL{E)G''{E)Tn{E)] (2) 

as G = GoT{Ep). This formula was originally derived 
for non-interacting electrons, but is in fact valid for inter- 
acting electrons in the low-bias limits. We have verified 
that this is indeed fulfilled to high numerical accuracy in 
our GW calculations by comparing the conductance ob- 
tained from Eq. ([1]) for small finite voltages with T{Ef) 
evaluated in equilibrium. 

The retarded Green function of the extended molecule 
is calculated from 

G^{E) = [{E + i7j)S-Ho + - ^Vii[G] 

-EL(^)-Sfl(£;)-S.c[G](S)]~i (3) 

where 7y is a positive infinitesimal and 

S.J = {(I^A^j) (4) 
Ho,tj = - + i>ion(r) + VH{r) + v^c{r)\(l>j}{5) 

Kc,y = (0>:Ec(r)|0j) (6) 

denote the overlap matrix, Kohn-Sham Hamiltonian and 
exchange-correlation potential, respectively. The matri- 
ces are evaluated in terms of a basis consisting of nu- 
merical atomic orbitals^, and are obtained from a DFT 
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supercell calculation performed with the real-space pro- 
jector augmented wave method GPAW^^. The electrode 
self-energies Y^LjR are obtained from the Kohn-Sham 
Hamiltonian of a bulk DFT calculation using standard 
techniques''^. 

The boundary conditions in the plane normal to the 
transport direction enter only via the electrode self- 
energies which are constructed from the electrode surface 
Green function^. The latter should represent an infinite 
surface but is here approximated by that of a periodic 
supercell with 4x4 Au atoms in the surface plane. We 
have found that this is a very good approximation when 
the surface Green function of the 4x4 cell is evaluated at 
a general fc_L-point (we use fcj^ = (0.1, 0.2) in coordinates 
of the surface Brillouin zone basis vectors). Using high 
symmetry points, in particular the Gamma point, can be 
problematic^. 

The term AT^J is the deviation of the Hartree potential 
from the groundstate DFT Hartree potential contained 
in i?o, 

AyH,.j = 2 ^ v,,M{eu\G\ ~ qI,\G^\) (7) 

kl 

In this expression q = —i J G^{E)dE and g'^ = 
—i J GQ{E)dE are the interacting and Kohn-Sham den- 
sity matrices, respectively, and 

is the bare Coulomb interaction in the atomic orbital 
basis. The factor 2 is due to spin. Ref. IH describes how 
the all electron Coulomb elements are obtained within 
the PAW formalism. The last term in Eq. ([3]) is the 
many-body exchange-correlation self-energy which in this 
work can be either the bare exchange potential or the GW 
self-energy. We note that setting AVh = and S^c = 14c 
in Eq. ([3]) we obtain the Kohn-Sham Green function, Gq. 



A. GW self-energy 

The (time-ordered) GW self-energy is given by 

SgW,^,(t) =iY,GkliT+)Wk^,lJiT) (9) 
kl 

where W is the screened interaction and the indices 
i,j,k,l run over the atomic basis functions of the ex- 
tended molecule. The screened interaction is calculated 
in the frequency domain as the matrix product W{uj) = 
e~^{uj)v with the dielectric function, e{uj) = 1 — vP{u)), 
evaluated in the random phase approximation. The irre- 
ducible density reponse function is computed in the time 
domain, 

P^0M{^) = -2iG,k{T)Gl,{^T). (10) 



where the factor 2 accounts for spin. Setting P = 
yields the Hartree-Fock approximation which thus corre- 
sponds to complete neglect of screening or equivalently 
complete neglect of correlations. Note that Eqs. ([ MTU)) 
involve time-ordered quantities defined on the Keldysh 
time contour. For completeness we provide the expres- 
sions for the real time components in Appendix |X] and 
refer the reader to Ref. 35 for more details. 




FIG. 1: (Color online) Schematic of a molecular junction with 
the "molecule" and "extended molecule" regions indicated 
(small and large box, respectively). The GW self-energy (E) 
is calculated on the molecule while the polarization (P) and 
screened interaction {W) are evaluated for the extended re- 
gion to ensure proper treatment of non-local screening effects. 
The electrodes and electrode-molecule coupling is described 
at the DFT level. A second order Feynman diagram is shown. 

As explained above, the matrices P, W , and Sow are 
calculated for the extended molecule. On the other hand 
it is clear that VF, and thus SgWi should have contribu- 
tions from polarization diagrams outside this region, see 
Fig. [1] Physically these diagrams describe the poten- 
tial acting on an electron propagating on the molecule 
due to the polarization that it induces in the lead. For 
this reason the self-energy will not be fully converged at 
the ends of the extended molecule region. To overcome 
this problem we only use the part of Sqw correspond- 
ing to the molecule and replace the remaining parts of 
the xc self-energy of the extended molecule by the DFT 
xc-potential. Symbolically, 

V^c ^GW{E) V^c , (11) 
'^xc ^xc ^xc J 

We stress that although we only include Sqw on the 
molecule, the interactions between electrons on the 
molecule and electrons in the electrodes (leading e.g. to 
image charge renormalization of the molecular levels) are 
included via diagrams of the form shown in Fig. [TJ We 
also note that the form ([TT|) implies that all metal atoms, 
both those inside and outside the extended molecule, are 
consistently described at the same (DFT) level. For the 
size of the extended molecule we have found it sufficient 
to include the gold atoms which are nearest neighbors to 
the sulfur or nitrogen atoms, i.e two gold atoms for the 
tip structures and six for the fiat structure depicted in 
Fig. [21 see Appendix [BJ We expect that this rather local 
screening response is special for covalent metal-molecule 
bonds. 
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B. Time/frequency dependence 

The time/frequency dependence of G, P, W, and Sqw 
is represented on a uniform grid ranging from —200 to 200 
eV with a grid spacing of Aw = O.OleV. We have veri- 
fied that the results are coverged with respect to both 
the size and spacing of this grid. Fast Fourier trans- 
form is used to switch between energy and time domains 
to avoid convolutions. The calculations are parallelized 
both over basis functions and over the time/frequency 
grid points. One should always have 77 > Auj to ensure 
proper representation of possible bound states. However, 
we have found that the conductance, and more generally 
the transmission function at any given energy, can be 
linearly extrapolated to the 77 = 0^ limit. This extrapo- 
lation has been performed for all the results presented in 
this work. 

The memory requirements for the GW calculations 
(defined mainly by the size of the P and W matrices) 
are approximately a factor 50 larger than for a corre- 
sponding DFT calculation. The GW calculations for the 
benzene junctions considered in the present work were 
performed in parallel on 100-400 cores and took about 
2 hours per selfconsisistency iteration. In comparison a 
DFT calculation for the same system took around 5 hours 
on 8 cores. 



C. Product basis 

The calculation of all of the Coulomb matrix elements, 
Vij^ki, is prohibitively costly for larger basis sets. For- 
tunately the matrix is to a large degree dominated by 
negligible elements. To systematically define the most 
significant Coulomb elements, we use the product basis 
technique of Aryasetiawan and GunnarssonS^. In this 
approach, the pair orbital overlap matrix 

Sij,ki = {nij\nki), (12) 

where nij{r) — (p* {r)ipj{r) is used to screen for the sig- 
nificant elements of v. 

The eigenvectors of the overlap matrix Eq. rep- 
resents a set of "optimized pair orbitals" and the eigen- 
values their norm. Optimized pair orbitals with insignif- 
icant norm must also yield a reduced contribution to the 
Coulomb matrix, and are omitted in the calculation of 
V. We have found that the basis for v can be limited to 
optimized pair orbitals with a norm larger than IO^^Oq 
without sacrificing accuracy. This gives a significant re- 
duction in the number of Coulomb elements that needs 
to be evaluated, and it reduces the matrix size of P{uj) 
and W{uj) correspondingly, see Appendix El 

D. Valence-core exchange 

Since both core and all-electron valence states are 
available in the PAW method, we can evaluate the contri- 



bution to the valence exchange self-energy coming from 
the core electrons. As the density matrix is simply the 
identity matrix in the subspace of atomic core states, this 
valence-core exchange reads 

core 
n 

where i,j represent valence basis functions and n rep- 
resent atomic core states. This contribution is added 
to Eqw describing the valence-valence interactions. We 
limit the inclusion of valence-core interactions to the ex- 
change potential, neglecting it in the correlation. This is 
reasonable, because the polarization bubble, P, involv- 
ing core and valence states will be small due to the large 
energy difference and small spatial overlap of the valence 
and core states. In general we have foimd that the effect 
of on molecular energy levels can be up to 1 eV— . 
For the benzene-like molecules considered in this work 
the effect is generally less than 0.4 eV for the frontier 
orbitals. 



E. Self-consistency 

Since Sqw a-nd AVu depend on G, the Dyson equation 
([5]) must be solved self-consistently in conjunction with 
the self-energies. In practice, this self-consistency prob- 
lem is solved by iteration. We have found that a linear 
mixing of the Green functions, 

GUE) = (1 - a)Gl-\E) + aG':^,{E) (14) 

with a = 0.15 generally leads to selfconsistency within 
20-30 iterations. 

Fully selfconsistent GW calculations are not standard, 
and in fact only few previous calculations of this type 
have been reported^^i^i^. Conventional GW bandstruc- 
ture calculations typically apply a one-shot technique 
where the self-energy is evaluated with a non-interacting 
Green function, Go, usually obtained from DFT~. In 
comparison, the selfconsistent approach has the imme- 
diate advantage of removing the Go-dependence, i.e. it 
leads to a unique solution. 

For an approach, like the present, where the chemical 
potential is fixed by the external boundary conditions, 
some kind of selfconsistency (though not necessarily the 
full GW selfconsistency employed here) is essential to en- 
sure charge neutrality. This is particularly important 
for cases where a molecular resonance lies close to the 
chemical potential. A shift in the energy of such a reso- 
nance could lead to a large change in its occupation. In a 
selfconsistent calculation this shift would be counterbal- 
anced, mainly by a change in the Hartree potential. On 
the other hand, the one-shot approach does not account 
for this effect and unphysical level alignments could occur 
as a result. 

Finally, the fully selfconsistent GW approximation is a 
conserving approximation in the sense of Baym^^. This 



5 



becomes particularly important in the context of trans- 
port where it ensures that the continuity equation is 
satisfied^^i^. We mention that the recently introduced 
quasi-selfconsistent GW method (not to be confused with 
the fully selfconsistent GW approximation used here), 
in which Go is chosen such as to mimick the interact- 
ing Green function as closely as possible, have shown 
that selfconsistency in general improves the band gaps 
of semiconductors as compared to standard one-shot 
calculations.— 



III. RESULTS 

In this section we discuss the results of our self- 
consistent GW calculations for the electronic struc- 
ture and conductance of the prototypical gold/ 1,4- 
benzenedithiolate (BDT), -benzenedithiol (BDT-^H), 
and -benzenediamine (BDA) molecular junctions. We 
argue that the thiol structure can be considered as a 
simple model for more complex gold/thiolate structures 
which have been proposed recently^^ but which are cur- 
rently too large to be treated satisfactorily at the GW 
level. The transport results are rationalized by consider- 
ing the alignment of molecular energy levels in the junc- 
tion. Here we show that both DFT and Hartree-Fock 
provide quantitatively and qualitatively wrong results by 
predicting a gap opening rather than reduction when the 
molecule is attached to electrodes. Finally, we investigate 
to what extent the GW results can be reproduced by a 
simple scissors operator applied to the DFT Hamiltonian. 



A. Junction geometries 

The junction geometries were optimized using the real 
space projector augmented wave method GPAW— . We 
used a grid spacing of 0.2 A and the PBE functional for 
exchange and correlation^^. The molecules were attached 
to Au(lll) surfaces, modelled by a seven layer thick 4x4 
slab, either directly (in the case of BDT) or via tips (in 
the case of BDT+H and BDA) as shown in Fig. H The 
surface Brillouin zone was sampled on a 4 x 4 /c-point 
grid, and the structures including molecule, Au tips, and 
outermost Au surface layers were relaxed until the resid- 
ual force was below 0.05 eV/A. The three structures are 
shown in the upper panel of Fig. [2] and some key bond 
lengths are given here^S. 

It is generally accepted that the hydrogen atoms dis- 
sociate from the thiol end groups forming a gold-thiolate 
structure^. Nevertheless, our total energy calculations 
show that the benzenedithiol structure has a slightly 
lower energy than the benzenedithiolate when inserted 
between two gold tips as shown in Fig. ^middle). In 
these calculations the hydrogens were either taken to be 
in the gas-phase or are adsorbed on the Au tips. In both 
cases the energy gain is less than 0.1 eV at the equi- 
librium distance but increases to 0.3 eV for a junction 



stretched by 1 A. We stress that our calculations do not 
include effects of entropy which becomes relevant at finite 
temperature, and furthermore assumes that hydrogen in 
the gas phase is the proper reference for the solvated 
proton and an electron in the electrode, i.e. the reaction 
H^+c^ -H- H2(g) is in equilibrium (in electrochemistry 
language we assume the standard hydrogen potential at 
pH=0). For these reasons our calculations are not suffi- 
cient to address the relative stability of thiols vs. thio- 
lates under the relevant experimental conditions. 

Importantly we note that new experimental evidence 
for the chemical structure of the gold-thiolate inter- 
face at the Au(lll) surfaceiS^J^ or at Au nanocluster 
surfaceai^ii^ has recently emerged, pointing to the ex- 
istence of polymeric SR-Au(I)-SR units where the for- 
mally oxidised Au(I) adatoms are chemically bound to 
thiolates and form a part of a more complex structure 
(see Ref. U3 and references therein). These complexes 
are currently too challenging to treat satisfactorily at the 
GW level. However, we have found that the electronic 
structure of such complexes is quite similar to that of the 
dithiol structure, see Sec. IIII Gl This is because the hy- 
drogen atoms play a role similar to that of the Au adatom 
in passivating the sulfur atoms. Therefore the transport 
properties of the dithiol structure should also be similar 
to those of the SR-Au(I)-SR units. 



B. Energy levels of isolated molecules 

A natural requirement for a method intended to de- 
scribe the energy levels of molecules in contact with elec- 
trodes, is that it should be able to describe the energy 
levels of isolated molecules. As we show below, the DFT 
approach fails completely in this respect underestimating 
the HOMO-LUMO gaps of the three molecules by 5-6 eV 
while GW energies lie within 0.5 eV of the target values. 

The gas-phase molecular structures have been relaxed 
in a 16 A cubic cell using GPAW grid calculations as 
described in the previous section. For consistency, all en- 
ergy levels have been calculated using the same double- 
zeta (DZ) atomic orbitals basis set. This is the same 
basis as used for the molecules in the transport calcula- 
tions presented in the next section. For the DFT and 
Hartree-Fock calculations we have found that the ener- 
gies obtained with the DZ basis agree with accurate grid 
calculations to within 0.2 eV. For the GW calculations 
the energies are within 0.1 eV of those obtained with a 
DZ-|-polarization basis set. 

Table U summarizes the results for the HOMO and 
LUMO orbital energies obtained from the DFT-PBE 
eigenvalues, selfconsistent Hartree-Fock, and selfconsis- 
tent GW. The energy levels have been identified as the 
peaks in the spectral function Tr[ImG""(£^)] extrapolated 
to 77 = 0"*" . Due to lack of accurate experimental data we 
have also performed DFT-PBE total energy calculations 
for the neutral, cation, and anion species to obtain the 
addition/removal energies (last column). This approach 
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FIG. 2: (Color online) Atomic structure of the BDA (left), BDT+H (middle), and BDT (right) molecular junctions. The 
lower panels show the transmission functions calculated using DFT-PBE (red) , Hartree-Fock (blue) , and the selfconsistent GW 
approximation (green). The insets show a zoom of the transmission functions around the Fermi energy (set to zero). The grey 
boxes indicate the conductance windows 0.007Go ± 50% and O.OlGo ± 50% which cover the the experimental values reported 
in Refs. and [^-d) respectively. 



has been shown to produce very accurate estimates of 
the experimental ionization and affinity levels of small 
molecules^. 

Relative to this reference, the DFT eigenvalues under- 
estimate the HOMO-LUMO gap of all three molecules 
by 5 — 6 eV, Hartree-Fock overestimates it by 2 — 3 eV, 
while the gap obtained with GW lies within 0.3 eV. These 
trends are consistent with a recent study of ionization 
potentials of a large number of molecules^. The main 
reason for the large underestimation of the gap by DFT is 
the presence of self-interactions in the FEE functionali^ 
On the other hand Hartree-Fock is self-interaction free; 
here, by virtue of Koopmans' theorem, the overstimation 
of the gap can be seen as a result of neglect of orbital re- 
laxations. The effect of the latter is included in GW via 
the screened interaction and this reduces the gap relative 
the unscreened Hartree-Fock result. 

We furthermore note that DFT places the LUMO of 
BDA and BDT-I-H below the vacuum level thus incor- 
rectly predicting the anions to be stable. For BDT, the 
LUMO level is predicted to be negative by all methods 
indicating the radical nature of this species. We note 
that our GW energies for BDT are in good agreement 
with previous MP2 calculations^. 



C. Conductance calculations 

The transmission functions of the relaxed junction ge- 
ometries were calculated as described in Sec. [TT] us- 
ing three different approximations for E^jc, namely the 
PBE xc-potential, the bare exchange potential (leading 
to Hartree-Fock theory), and GW. The former choice 
corresponds to the standard DFT-approach. All calcu- 
lations employ a double-zeta basis set for the molecules 



Molecule 


Orbital 


DFT-PBE 


HF 


GW 




BDA 


HOMO 


-4.1 


-7.2 


-6.2 


-6.8 


(CsHgNa) 


LUMO 


-0.9 


3.9 


2.9 


2.3 




H-L Gap 


3.2 


11.1 


9.1 


9.1 


BDT+H 


HOMO 


-5.1 


-8.0 


-6.9 


-7.5 


(C6H6S2) 


LUMO 


-1.3 


3.3 


2.2 


1.3 




H-L Gap 


3.8 


11.3 


9.1 


8.8 


BDT 


HOMO 


-5.7 


-8.6 


-7.9 


-8.3 


(C6H4S2) 


LUMO 


-5.1 


-1.6 


-2.3 


-2.7 




H-L Gap 


0.6 


7.0 


5.6 


5.6 



TABLE L Calculated frontier orbital energies of the molecules 
in the gas-phase. All energies are in eV and measured relative 
to the vacuum level. DFT-PBE refers to the Kohn-Sham 
eigenvalues while ABtot represents addition/removal energies 
obtained from self-consistent total energy calculations of the 
neutral, anion and cations at the DFT-PBE level. 



and double-zeta with polarization for the Au atoms. The 
results are shown in the lower panels of Fig. [2l and the 
corresponding conductances are summarized in Table [Hi 

The conductance of BDA and BDT-f-H calculated with 
the fully selfconsistent GW approximation agree well 
with the experimental values reported in Refs. yand^for 
benzenediamine and Refs. y-ly for benzenedithiol as in- 
dicated by the grey boxes in Fig. [^l^left -I- middle). In con- 
trast, DFT and Hartree-Fock respectively overestimates 
and underestimates the experimental conductances by 
factors 5-10. Our DFT result for the BDA junction is 
in good agreement with previous calculations^"^'^^. 

It is striking that the conductance of the "classical" 
BDT junction shown in Fig. [21[right), is predicted by all 
three methods to be significantly higher than the experi- 
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Method 


BDA 




BDT+H 


BDT 


DFT-PBE 


2.1 ■ 10' 


2 


5.4 ■ 10"^ 


2.8 ■ 10"^ 


HF 


4.0 • 10" 


4 


2.7- 10"^ 


5.7- 10-^ 


GW 


4.2 • 10" 


3 


1.0 ■ 10"^ 


8.3 • 10"^ 


GoWo(PBE) 


8.0 • 10" 


3 


1.6 • 10"^ 


7.5 • 10"^ 


GoWo(HF) 


2.2 • 10" 


3 


9.8 • 10"^ 


8.7-10"^ 



TABLE II: Calculated conductance in units of Go for the 
three junctions shown in Fig. (2] The last two rows refer to 
non-selfconsistent GW calculations based on the DFT-PBE 
or Hartree-Fock Green function, respectively. 

mental value (we obtain the same high conductances for 
BDT between tips, i.e. the structure in Fig. [JJmiddle) 
without hydrogen on sulfur). Our DFT conductance is in 
overall good agreement with the large number of previous 
calculations for the same or similar similar structure^^. 
The high conductance is clearly due to a strong peak in 
the transmission function close to the Fermi level. The 
peak moves to higher energies when going from DFT- 
PBE over GW to Hartree-Fock, opposite to the trend 
normally seen for occupied states. The peak comes from 
an unsaturated p-orbital on the sulphur atoms and is dis- 
cussed in more detail in Sec. IIII Gl 

These results suggest that the structures probed in ex- 
periments on benzenedithiol junctions involve a chemi- 
cally passivated form of the thiolate linker group. As 
we show in Sec. IIII G[ the high conductance of BDT 
is due to unsaturated p-states on the sulphurs with en- 
ergy close to the Fermi level (the energy of this orbital 
does not change considerable with DFT, Hartree-Fock 
and GW). In the thiol and SR-Au(I)-SR structures, these 
states form bonds to H and the Au adatom, respectively, 
and are thereby shifted away from the Fermi level. On 
the other hand, the electronic structure and transmis- 
sion functions of BDT+H and SR-Au(I)-SR structure 
are rather similar indicating that the BDT-I-H can be 
viewed as a simple model of the more complex SR-Au(I)- 
SR structure. 

It should of course be kept in mind that experiments 
are performed in solution and at room temperature and 
are subject to variations in the detailed atomic struc- 
ture. Thus the measured conductance values should not 
be considered as highly accurate references for theoretical 
calculations on idealized junctions. 

D. Energy level alignment 

In Fig. |3]we show the calculated HOMO and LUMO 
energy levels of BDT-I-H and BDA in the gas-phase and 
in the junction. All energies have been aligned relative to 
the vacuum level. At this point we note that an accurate 
description of the vacuum level, i.e. the work function, 
can in general be difficult to obtain with an atomic orbital 
basis^i. However, by using more diffuse basis functions 
(an energy shift of 0.01 eV has been used for all Au basis 




FIG. 3: (Color online) Upper panel shows the change in the 
HOMO-LUMO gap as the molecules are brought from the 
gas-phase into the junction. Lower panel shows the energy of 
the HOMO and LUMO levels (relative to vacuum) in the gas 
phase and in the junction. The left and right set of levels cor- 
respond to BDT+H and BDA, respectively. The calculated 
work function of the gold junction (with tips) in absence of a 
molecule is $ = 4.9 eV as indicated. 



functions throughout this work^^) we obtain a work func- 
tion for Au(lll) of 5.4 eV in good agreement with the 
experimental value of 5.31 eV— . At the position of the 
molecule, i.e. in the region between the two tips, the elec- 
trostatic potential from a calculation where the molecule 
has been removed, converges to a constant value of 4.9 eV 
above the metal Fermi level. This value has been used as 
reference for the vacuum level in Fig. |31 In the junction 
where the levels are broadened due to hybridization with 
the metal states, the level positions have been defined as 
the first moment of the projected density of states of the 
relevant molecular orbital, Im(V'n|G'"(i?)|V'n)- Here the 
\ipn) are obtained by diagonalizing the KS Hamiltonian 
within the molecular subspace. 

The orbital energies obtained from a GW calculation 
include the dynamical response of the electron system 
to the added electron/hole via the correlation part of 
the self-energy. In general the correlations will shift the 
occupied levels up and the empty levels down relative 
to the bare Hartree-Fock energies. When a molecule is 
brought from the gas-phase into a junction the electronic 
character of its environment changes from insulating to 
metallic. The enhanced screening should thus cause the 
gap to shrink (neglecting shifts due to hybridization) as 
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compared to its gas-phase value. However, it has recently 
been shown for molecules weakly bonded to a metal sur- 
face, that this effect is completely missing in effective 
single-particle theories based on a (semi)local description 
of correlations^-SiS. 

As a result, in our DFT and Hartree-Fock calculations, 
the change in the frontier orbital energies induced by the 
coupling to gold is completely governed by the effect of 
hybridization which tends to open the gap by 0.7 eV 
for both molecules, see top panel of Fig. |3l In contrast, 
the GW gap is reduced by 1 eV due to the enhanced 
screening in the contact. Since the hybridization shift is 
of course also present in he GW calculation, we conclude 
that the enhanced screening due to the metal contacts 
reduces the HOMO-LUMO gap by 1.7 eV relative to 
the value in the gas-phase. 

Note that we refrain from using the term "image 
charge effect" to describe the renormalization of molec- 
ular orbitals. This term is appropriate for weakly cou- 
pled molecules where the screening of the added elec- 
tron/hole takes place within the metal. For intermediate 
or strongly coupled molecules, there is no clear distinc- 
tion between metal states and molecular states, and the 
screening is more appropriately described as dynamical 
charge transfer-^. 

From Fig. H it follows that the HOMO level of the 
molecules in the junction is predicted by DFT-PBE to 
lie 0.5 - 0.7 eV higher than obtained with GW. This 
agrees well with a recent study of benzenediamine deriva- 
tives on gold(lll) which showed that DFT places the 
HOMO level about 1 eV too high compared to UPS 
measurements^^. The fact that the DFT-PBE descrip- 
tion of the energy levels is better for the adsorbed 
rather than isolated molecules may be seen as a result of 
the metallic screening build into the DFT xc-functional 
via its origin in the homogeneous electron gas^l. It 
should also be noted that the error (compared to GW) 
of the DFT eigenvalues is significantly larger for the 
LUMO than for the HOMO. This is in good agreement 
with previous plane wave calculations for molecule/metal 
interfaces^. 



E. Go Wo calculations 

To test the role of selfconsistency (in the GW 
and Hartree self-energies) we have performed non- 
selfconsistent GqWo calculations using both the DFT- 
PBE Green function and the (self consistent) Hartree- 
Fock Green function as the initial Gq. The results for 
the conductance of all systems are summarized in the last 
two rows of Table [Hi We notice is that the conductance 
can vary by a factor of three depending on Gq. We also 
note that the Hartree-Fock starting point comes closer 
to the selfconsistent result. This is because the Hartree- 
Fock Green function is an overall better approximation 
to the GW Green function than is the DFT Green func- 
tion (see e.g. the comparison of frontier orbital energies 



in Fig. |3]). As an example Fig. |4] shows the calculated 
transmission functions for the BDA junction in a larger 
energy range around Ep. 

For the BDA and BDT+H junctions, GoWo[PBE] 
overestimates the conductance while GoWo[HF] under- 
estimates the conductance relative to GW. This can be 
understood as follows: Since DFT-PBE (Hartree-Fock) 
underestimates (overestimates) the HOMO-LUMO gap, 
the use of these Green functions to evaluate the GW self- 
energy will lead to an overestimation (underestimation) 
of the screening. As discussed above the screening con- 
tained in the correlation part of the GW self-energy tend 
to reduce the HOMO-LUMO gap. This reduction of the 
gap is thus overestimated in the Gq Wq [PBE] calculation 
and underestimated in the GoWo[HF] calculation which 




FIG. 4: (Color online) Transmission functions for the 
BDA junction calculated with selfconsistent GW and non- 
selfconsistent Go Wo using either the Hartree-Fock or the 
DFT-PBE Green functions as input. 

The Go Wo results for the strongly coupled BDT junc- 
tion show less variation. This is perhaps surprising given 
the large difference between the DFT-PBE and Hartree- 
Fock results shown in Fig. [51 However, DFT-PBE and 
Hartree-Fock give almost equal density if states at the 
Fermi level (the transmission functions at the Fermi en- 
ergy are also rather similar), and therefore the screening 
contribution obtained with the two choices for Go is also 
very similar. 

F. Scissors operator 

In this section we investigate to what extent the GW 
transmission function can be reproduced by a DFT cal- 
culation where the occupied and unoccupied molecular 
levels have been shifted rigidly to match the GW ener- 
gies. This is illustrated by applying a scissors operator 
(SO) to the Au-BDA-Au junction shown in Fig. [IJleft). 
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After shifting the molecular levels we perform a non- 
selfconsistent calculation of the transmission function. 
A similar procedure has previously been successfully 
used to include image charge effects and correct for self- 
interaction errors in DFT-transport calculations^i^iS. 
We refer the reader to Ref. [23| for more details on the 
SO technique applied here. 




Shift of HOMO (eV) 

FIG. 5: (Color online) DFT conductance of the Au-BDA-Au 
junction for different shifts of the occupied and unoccupied 
molecular orbital energies. The three iso-contours with val- 
ues of 4.2x10-3, 6.4x10-3 and 2.1x10-2 in units of Go corre- 
spond to the GW, experimental and DFT conductance values, 
respectively. The dot indicates the conductance obtained by 
fitting the DFT HOMO and LUMO level positions to the cor- 
responding GW level positions. Inset: Transmission function 
calculated with GW (green) and DFT-hSO (red). The SO 
shifts are -0.6 eV and 3.8 eV for the occupied and unoccupied 
molecular orbital, respectively, which makes the HOMO and 
LUMO energies of the DFT calculation coincide with the GW 
energies. 

In Fig. [S]we show contour plots of the DFT conduc- 
tance for the BDA junction where the shift of the occu- 
pied and unoccupied molecular orbital energies has been 
varied independently over 4 eV. The three values for the 
contour lines shown correspond to the conductance ob- 
tained with GW (0.0042Go), the experimental conduc- 
tance (0.0064Go) and the DFT conductance (0.021Go), 
respectively. The dot indicates the energy shifts which 
make the HOMO and LUMO levels of the DFT calcu- 
lation match the corresponding levels of the GW cal- 
culation; the required SO shifts are -0.6eV and 3.8eV 
for the occupied and unoccupied molecular orbitals, re- 
spectively. Shifting the levels by this amount reduces 
the DFT conductance by a factor of 3 from 0.02Go to 
0.0072Go. Interestingly, the GW conductance is not re- 
produced by these shifts; it is even lower by a factor of 
1.5. In fact, to reproduce the GW conductance a shift 
of about -l.SeV of the DFT HOMO is required (keep- 
ing the LUMO position fixed at the GW position). This 
shows that while the renormalization of the molecular 



level energies can explain the main part of the difference 
between the DFT and GW conductance, the different 
shape of the transmission resonances also plays a role. 
This is clear from the inset which shows the GW trans- 
mission function (green) and the DFT transmission with 
SO chosen to match the GW HOMO and LUMO levels 
(the dot in the main figure). The lower conductance ob- 
tained with GW is seen to be a consequence of a faster 
decay of the HOMO resonance towards the Fermi level. 
This difference comes from the energy dependence of the 
GW self-energy. 



G. Thiol vs. thiolate structures 

In Fig. ini we compare the DFT transmission func- 
tions for: (a) the "classical" structure of benzenedithi- 
olate [structure in Fig. [^fright)] (b) benzenedithiol be- 
tween tips on Au(lll) [structure in Fig. OJmiddle)], and 
(c) benzenedithiolate in a SR-Au(I)-SR molecular unit 
form [structure 2 of Fig. 1 in Ref. |T5[. The conductance 
(essentially the transmission function at the Fermi en- 
ergy) obtained for structures (b) and (c) is very similar 
but markedly different from (a). The high conductance 
of the structure (a) is due to the strong transmission res- 
onance lying just below the Fermi level. This resonance 
has also been found in many previous studies and seems 
to be a characteristic and robust feature of this junc- 
tion. As shown in this work this transmission resonance 
is also present in GW and Hartree-Fock calculations. In 
contrast, for both structure (b) and (c) the transmission 
function is rather flat in an energy window of ±1 eV 
around the Fermi level. 

To examine the origin of the different transmission 
functions found for structures (a) and (b,c) we found it 
useful to consider the molecular orbitals of the C6H4 moi- 
ety of the three molecules. By diagonalizing the (Kohn- 
Sham) Hamiltonian of this part of the molecules we found 
that the orbital depicted to the left in the upper panel of 
Fig. El and which constitute the HOMO-1 of the C6H4 
moiety, is responsible for all the structure in the transmis- 
sion functions below the Fermi level. Note, that the or- 
bitals obtained in this way are different from the HOMO 
and LUMO levels shown in Fig. [3] which were obtained 
by diagonalizing the Hamiltonian for the entire molecules 
including the SH and NH2 end groups. 

The different panels of Fig. [7] show the projected den- 
sity of states (PDOS) of the HOMO-1 for the three struc- 
tures together with the PDOS of the sulfur p-orbital 
to which the HOMO-1 couples. Within the Newns- 
Anderson model2^, the sulphur p-orbital is the so-called 
group orbital. Note, that the PDOS of the p-orbital 
has been calculated in the absence of coupling to the 
HOMO-1 as it should in the Newns-Anderson model. 
The transmission function is then essentially the prod- 
uct of the PDOS of the HOMO-1 and the PDOS of the 
group orbital2S. It is clear that the origin of the (double) 
transmission peak around -1 eV for structure (a) is due 
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FIG. 6: (Color online) DFT-PBE transmission functions for 
three different structures of the Au-BDT-Au junction: a) the 
"classical" benzenedithiolate, b) benzenedithiol, and c) the 
SR-Au(I)-SR complex (with two benzene molecules replaced 
by CH3 units for simplicity). The transmission function of b) 
and c) are rather similar while that of a) shows a peak close 
to the Fermi level. 



to a resonance formed by the HOMO-1 and the sulfur 
p-orbital. The chemical passivation of sulfur, by either 
hydrogen or the Au adatom, implies that the PDOS of 
the p-orbital splits into bonding and anti-bonding states 
around -3 eV and 3 eV, respectively. This in turn shifts 
the PDOS of the HOMO-1 down in energy. In particu- 
lar its magnitude around the Fermi level is lowered and 
as a consequence the transmission function (being essen- 
tially the product of the two curves) is suppressed in a 
broad energy window around Ep. Thus chemical pas- 
sivation of the sulfur is the main reason for the lower 
conductance observed in structures (b,c) as compared to 
(a). A secondary effect, giving rise to differences in the 
transmission of structure (b) and (c) is the different in- 
terface dipoles which shift the electrostatic potential at 
the C6H4 moiety by different amounts. This shift, how- 
ever, has little influence on the conductance due to the 
flatness of the transmission functions around Ep. 

To verify this scenario, we have applied a scissors op- 
erator of 1 eV to the C6H4 moiety of structure (b) in 
order to align the onsite energy of the HOMO-l to that 
of structure (c). The resulting transmission function (not 
shown) is very similar to that of (c) and the conductance 
is 3.5 • 10~^Go compared to 3.3 • 10~^G'o obtained for 
structure (c). On the other hand, the conductance of 
BDT cannot be brought below O.IGq by sifting the levels 
of the C6H4 moiety down (by up to 2 eV) because of the 



strong peak around -1 eV due to the sulphur p-orbital. 

Although the above picture is based on the DFT elec- 
tronic structure, the qualitative similarities of the DFT 
and GW transmission functions in Fig. [2] implies that 
the same picture applies to the GW calculations. 
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FIG. 7: (Color online) Projected density of states for the 
HOMO-l of the C6II4 moeity and its group orbital (sulfur p 
orbital). These orbitals essentially determine the transmission 
function around the Fermi energy for all the three structuers. 
It should be noticed that for structure (a) (the "classical" 
BDT structure) the sulfur p-orbital has a peak in the PDOS 
just below the Fermi level which is responsible for the high 
conductance. In (b) and (c) the sulfur is chemically passivated 
and the PDOS of the p-orbital splits into bonding and anti- 
bonding states at ±3 eV thereby lowering the conductance. 



IV. SUMMARY 

We have presented a first-principles method for mod- 
elling quantum transport in molecular nanostructures be- 
yond the single-particle approximation. The method is 
based on non-equilibrium Green functions and applies to 
the general case of a finite bias voltage, but in this work 
we focused on the zero bias regime. The conductance of 
benzenedithiolate (BDT), benzenedithiol (BDT-I-H), and 
benzenediamine (BDA) was calculated using the selfcon- 
sistent GW approximation. In contrast to standard DFT 
and Hartree-Fock methods, the GW approximation was 
found to yield consistently accurate values for the en- 
ergy levels of both isolated and contacted molecules due 
to its proper treatment of self-interaction and dynami- 
cal screening. In general, results obtained with GW for 
the electronic conductance and energy gaps of contacted 
molecules lie in between the values obtained with DFT 
and Hartree-Fock. The latter methods respectively over- 
estimates and underestimates the screening and none of 
them are able to detect the change in the molecule's elec- 
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tronic environment when it is coupled to electrodes. 

Non-selfconsistent GqWo calculations were found to 
yield conductance values within 50% of the GW results 
depending on the initial Gg. It was shown that the main 
difference between the GW and DFT calculations comes 
from the renormalization of the position of the molecular 
energy levels. However, changes in the shape of trans- 
mission resonances, and thus the conductance, also occur 
due to the energy-dependence of the GW self-energy. 

The GW conductance of benzenediamine was found 
to be in good agreement with experiments. For ben- 
zenedithiolate between Au(lll) surfaces we found a con- 
ductance well above O.IGq and thus conclude that this 
structure cannot be responsible for the measured conduc- 
tance around O.OIGq. On the other hand, a conductance 
close to this value was found for a hydrogenated ben- 
zenedithiol junction which, as we demonstrated, repre- 
sents a reasonable model of more complex gold/thiolate 
structures where the chemical passivation of sulfur is pro- 
vided by a gold adatom rather than hydrogen. These 
gold/thiolate structures are presently too demanding for 
our selfconsistent GW method, however, our results sug- 
gest that such structures are responsible for the O.OIGq 
conductance reported experimentally. 

In conclusion, we showed that a consistent and quan- 
titatively accurate description of energy level alignment 
and charge transport in phase-coherent molecular con- 
ductors can be obtained from first-principles with the 
(computationally feasible) GW approximation. We be- 
lieve this development is important for increased the syn- 
ergy between theory and experiments in molecular elec- 
tronics which is essential for continued progress in the 
field. 
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Appendix A: The GW self-energy 

Let U denote the rotation matrix that diagonalizes the 
pair orbital overlap Sij^ki = {nij\nki), i.e. WSU = al. 
The columns of U are truncated to those which have 
corresponding eigenvalues cr, > lO^^^a^^. We then only 
calculate the reduced number of Coulomb elements 

Vqq' = {nq\ ,^^,: \ng-), (Al) 



where riqir) are the optimized pair orbitals 

nqi'^) ^'^nij{r)U.,j,q/y^, (A2) 

which are mutually orthonormal, i.e. (ng|rtg') = Sqqi. 

The determination of the GW self-energy proceeds by 
calculating first the lesser and greater components of the 
polarization matrix in the time domain 

P^UW^^^G<{t)G>it)\ (A3) 

(A4) 

where the factor 2 accounts for spin and we have used 
G>(-t) = -G>{t)*. This is then downfolded to the 
reduced representation 

Pf,' = E V^^Ul^^.PtuUkUq'^. (A5) 

The screened interaction can be determined from the 
lesser and greater polarization matrices, and the static 



interaction Vqqi , via the relations 

p^{t) = e{t){p>{t)^p<{t)), (A6) 

W''{uj) = [/ - vP''{lo)]-^v, (A7) 

W>{uj) = W {uj)P> {uj)W (A8) 

W<{lo)^W>{lo)-W{uj) + W'\uj)\ (A9) 



where all quantities are matrices in the optimized pair 
orbital basis and matrix multiplication is implied. We 
obtain the screened interaction in the original orbital ba- 
sis from 

Wtu (^) « E U^^.'t^qWf^, {u,)^Uti^^. , (AlO) 

qq' 

which is an approximation due to the truncation of the 
columns of U . Finally the GW self-energy can be deter- 
mined from 

siw..,w = *EGiw^t<w (All) 

kl 

^hwit) - 0{t) (I]>w(t) - SgwW) + Sm.- (A12) 
The exchange and Hartree potentials are given by 

^,.,j^iY,^^k,iiG<{t = 0) (A13) 

kl 

Vn.,, = -2^J2y^,,klGtl{t ^ 0) (A14) 

kl 

The latter equals the first term in Eq. ([7|). 

In equilibrium the lesser and greater Green functions 
are obtained from 

G<iE) = -fiE-f,)[G'-{E) - G^-{E)^] (A15) 
G>{E) = (1 - f{E - fi))[G^{E) - G^{E)^] (A16) 

where / is the Fermi-Dirac distribution function. For 
self-consistent calculations, equation (|A3l) - (IA16p are it- 
erated untill convergence in G. 



Appendix B: Size of the "extended molecule" 

As explained in Sec. Ill Al the GW self-energy of the 
molecule has contributions from screening diagrams ex- 
tending into the electrodes. This is included by calculat- 
ing P and W for the extended molecule which includes 
the Au atoms closest to the molecule. In Fig. |S]we show 
the dependence of the transmission function on the size of 
the extended region. The latter has been varied from zero 
Au atoms (extended molecule coincides with molecule), 
to the Au tip atoms, to the four-atom Au pyramids on 
each side of the molecule. From this we conclude that 
it is sufficient to include polarization diagrams for the 
Au tip atoms as was done for the calculations presented 
in the main text. We mention that, for simplicity, the 
calculations shown in Fig. [5] have been performed at the 
GoWo(PBE) level and using a minimal single-zeta basis 
set. However, the conclusions regarding the convergence 
of the GW self-energy with respect to the size of the ex- 
tended molecule, should be equally valid for the case of 
selfconsistent GW with the DZ/DZP basis set used for 
all calculations in the main text. 
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FIG. 8: (Color online) Transmission function for the BDA 
junction calculated at the GoWo(PBE) level and with a min- 
imal single-zeta basis set. The three curves correspond to 
different number of gold atoms in the extended molecule re- 
gion. 
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